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Abstract. Fritz Gesztesy's varied and prolific career has produced many 
transformational contributions to the spectral theory of one-dimensional Schro- 
dinger equations. He has often done this by revisiting the insights of great 
mathematical analysts of the past, connecting them in new ways, and rein- 
venting them in a thoroughly modern context. 

In this short note we recall and relate some classic transformations that 
figure among Fritz Gestesy's favorite tools of spectral theory, and indeed 
thereby make connections among some of his favorite scholars of the past, 
Bohl, Darboux, and Green. After doing this in the context of one-dimensional 
Schrodingcr equations on the line, we obtain some novel analogues for discrete 
one-dimensional Schrodinger equations. 

Dem einzigartigen Fritz gewidmet. 



1. Introduction 

In 1906 |3j, Bohl introduced a nonlinear transformation for solutions of Sturm- 
Liouville equations, which is an exact, albeit implicit, counterpart to the Liouville- 
Green aproximation |23) . Bohl used the transformation as a tool in oscillation 
theory, and this has continued to be the main use of the Bohl transformation in the 
hands of later authors. Notably, Rab |24j used the Bohl transformation to prove 
necessary and sufficient conditions for oscillation of solutions, and showed that it 
is an effective foundation for Sturm-Liouville oscillation theory. Willett's lecture in 
[29j provides a clear description in English of the Bohl transformation in oscillation 
theory, including the contributions of Rab, while Reid's monograph |27j compares 
and contrasts it with the Priifer transformation. See also [261 1131 1161 1171 118) . 

In [10] §4, Davies and Harrell introduced a non-oscillatory variant of the Bohl 
transformation to connect the notions of Liouville- Green approximation, Green 
functions, and Agmon metrics for exponential decay of solutions. Some spectral 
bounds were derived as consequences. This analysis was extended in a series of 
articles by Chernyavskaya and Shuster (e.g., [H [5j [8]), to address questions of 
solvability, regularity, estimates of Green functions, and asymptotics in Sturm- 
Liouville theory. 

In this note we begin with a largely expository treatment of the classic Bohl 
transformation, concentrating for simplicity on the situation where all coefficients 
are real and regular and the Sturm-Liouville equation is in the standard form of 
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the one-dimensional Schrodinger equation. Then in the last section we show how 
the technique can be adapted to the case of a discrete Schrodinger equation on the 
integers. 

2. The interplay of the Bohl and Green functions 

Let V be real-valued and continuous, and consider a solution basis for the 
Sturm-Liouville equation 

(2.1) - u" + V{x)u = 0; 

we may normalize the basis u\^{x) so that W[u\,U2\ '■= u\vl 2 — u^vl-y = 1. (There 
is no assumption of an eigenvalue 0. For our purposes a possible nonzero spectral 
parameter has simply been incorporated into V .) The Bohl transformation maps 
this solution basis onto a second solution basis with remarkable properties, some of 
which are collected in a nutshell version in this section. 

Definition 2.1. Given a solution basis {ui^(x)} of (|2.ip . chosen so that the 
Wronskian W[iti,it2] := u\v! 2 — u^u'i — 1, we define the diagonal junction by 

(2.2) Z[ui,u 2 ](x) := ( Ul (x)u2(x)) 1/2 . 

The Bohl transformation of {ui^(x)} is an equivalent solution basis of (|2.ip . defined 
in terms of {11^2(2;)} by 

(2.3) B:{u 1 (x),u 2 (x)}^{<f ) ± (x):=Z(x)exp(±J 

Remark 1. a) The choice of the complex phase of the square root in 
is unimportant, but should be continuous in x. For brevity we write Z(x) for 
Z[u\,U2\{x) when the dependence on ux,2 is clear. The reason for calling it the 
diagonal function is that, as will be seen below, 

Z 2 (x) =G (x,x), 

where G (x, x) is the diagonal of a certain Green function G (x, y) for (I2.1[) . Among 
the useful properties of the function Z is that it solves the diagonal differential 
equation 

(2.4) J[Z] := -Z" + V(x)Z - ^ = 0, 

[Toi rm fi8i [29]. 

b). In fact, with the oscillatory situation in mind Bohl originally wrote the solution 
basis in the Liouville-Green form 

1 



sin 



R{t)dt) 



and 



R \J J 

1 



cos 
R 



R{t)dt , 



which is equivalent to (|2.3p under the identification 2Z 2 — > i/R and some harmless 
linear combinations. 

c). We recall that Gesztesy and Simon [12] have made connections between the 
Krein spectral shift function, the related Xi function, and the diagonal of the Green 
function. 
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Calling upon |10j §4, we collect some facts, which are verifiable directly: 
Theorem 2.1. 

(1) If ui(x) and U2{x) are solutions to (|2.1[) such that W[u\,U2\ — 1, and 
Ui(x)u2(x) does not vanish on the interval (a, 6), then Z[u\, U2\{x) satisfies 
(ESI) on (a,b). 

(2) If Z is a nonvanishing solution of (12.4[) on (a,b), then <jy^{x) as defined 
in (|2.3[) provide a pair of independent solutions of (|2.1j) on (a, b). In 
particular, each (jy^ix) is a linear combination of {u 1.2} and vice versa. 

(3) If x> t < := max(a;, y), resp. min(x,y), then 

(2.5) G (x, y) := Z(x)Z(y) exp (- f > -^j-dt 



is a Green function for (|2.1[) . in the sense that 

(2-6) (-£2 + Go(x ' y) = 6{ - X ~ v) - 

As an integral kernel, Go defines the inverse of a particular realization of — ^3 + 
V , but not a priori one for which the domain of definition includes u\ or U2, because 
of a possible mismatch of boundary conditions at finite points. This issue is not 
important for questions of oscillation or asymptotic behavior at infinity, but another 
concern remains, namely the possibility that Z vanishes, which would invalidate 
the transformation. Because of this we recall that in the absence of imposed finite 
boundary conditions, complex solutions can always be used to prevent Z from 
vanishing: 



Lemma 2.2. Suppose that u is a solution of (|2.ip on a finite or infinite interval 
(a, b), and that at some xq € (a, b), Re(w(xo))Im(u(xo)) 7^ and u'{xq)/u{xq) R. 
Then u does not vanish on (a, 6). 



PROOF. Because V is real- valued, Reu and Im« each satisfy (12. 1|) . and it 
therefore suffices to show that they are independent. 
Letting a = u' (xq) / u(xg) , a calculation shows that. 

Reu'(xo) Imu(xo) 
— - — Re a — Ima — - , 

Keu{xa) Reu(xo) 

Im u'(x ) _ Re a + Ima Re«(ip) 
Iuiu(xq) Imjt(io) 



It follows that 



Reu'(xo) Imu'(xo) T flmtt(xo) Rew(xo) 
- -lm a 



Reu(xo) Im«(io) \Rew(xo) Imu(xo) 

and therefore 

W[Imu,Reu] = -Ima ((Reu(a;o)) 2 + (Imn(i„)) 2 ) ^ 



□ 



This standard lemma implies that given any two linearly independent solutions 
of (|2.1[) it is always possible to find a pair of complex- valued linearly independent 
combinations that are non- vanishing on (a, b). The Wronskian of the new pair may 
be set to 1 by multiplying one solution by an appropriate constant, justifying the 
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conclusions of Theorem l2.ll It is of some use to consider a particular Z determined 
as follows. 

By construction any solution u as set forth in Lemma l2.2l will be linearly inde- 
pendent of its complex conjugate. We may therefore choose a complex number a 
so that W[u, a 2 u) = 1, which ensures that 

(2.7) Z{x) = a\u{x)\. 

Indeed, arg(a) is restricted to the values &7r/4 for integer k, as can be seen for 
example from (12.41) . which implies that 



(2.8) _| ur + FM = _i_ eR . 

(In passing we note the implication that the expression on the left does not change 
sign.) In fact, there are only two truly distinct cases for arg(a), viz., and \, 
due to the simple scalings in (|2.4j) and Theorem 12.11 when Z is replaced by iZ. If 
a > 0, i.e., Z(x) > 0, then the solutions ^ do not change sign, which corresponds 
to the case of disconjugacy for the ODE (12.11) in the classical theory [201 126) . This 
situation was the focus of [10] , in which some spectral bounds were derived and it 
was argued that 1/2Z 2 defines an Agmon metric. 

Otherwise, it may be assumed without of loss of generality that arg(a) = j, 
for which solutions may oscillate, in that their arguments increase or decrease by 
nit for n > 1 as s -> oo. A central question of Sturmian theory is whether solutions 
oscillate infinitely often, or only finitely often. When the increase in the argument 
of a solution is infinite, the equation (I2.1[) is said to be oscillatory. An approach to 
oscillation theory, equivalent to that of I24L 129] but bringing out the role of Green 
functions, can be based on the following version of a result of Gagliardo, as cited 
in 



Corollary 2.3 (Cf. [29], Corollary 3.2.). Suppose that (|2.1|) holds on an 
infinite interval (a, oo), and let Gb(x, y) be the Green function constructed according 
to the prescription leading to (|2.7p . Then either Gs{x,y) G M for all x,y, in which 
case the solution basis (j)^ is nonoscillatory, or else: The phase of has only finite 
increase on (a,oo) iff 1/G(x,x) € L 1 (a,oo). 

As has been known since the work of Wigner and von Neumann and Wigner, 
it is possible for eigenvalues to be embedded in the continuous spectrum of Sturm- 
Liouville equations I22L 125] . This phenomenon requires oscillatory solutions to be 
square-integrable, and can thus be related to the Bohl transformation as follows: 

Corollary 2.4 (Rab). An oscillatory solution of (|2.ip . written in the form 

exists and is square integrable if and only if the nonlinear equation 
—w"(x) + V(x)w(x) = — - 



Aw 3 (x) 

has a square-integrable solution. 

We close this section with a Darboux-type factorization [9], 114] . the novel 
feature of which is the role played by the diagonal function and the Bohl solution 
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basis (|2.3|) . For any complex valued, nonvanishing function Z(x) € AC 1 [(a, &)], 
define 

+r i d Z' 1 
V [Z] : =^-Z T 2^ 

It is immediate to see that 

v ± [z]<j> ± =Q, 

where (fr^ are defined in terms of Z by (|2.3p . A further calculation reveals that 



provided that Z satisfies the diagonal differential equation (|2.4|) . An alternative 
way to express the two factorizations in (|2.9[) is that 



-^2 +V(»)= (X>|Z| J 
where * designates the formal adjoint operation. 

3. The discrete form of the Bohl J[Z] 

In this section we show that most of the transformations and relationships pre- 
sented in the first section have counterparts for discrete one-dimensional Schrodinger 
equations. (Part of the material in this section has appeared in a preprint [19) . 
which has been expanded and divided for publication as two articles.) some de- 
tails are rather different from the continuous case, making it uncertain how far 
the analogy goes, especially in the oscillatory case. A full-fledged oscillation the- 
ory for discrete problems based on an analogue of the Bohl transformation and its 
connection to Green functions would be an interesting next project. 

Let A denote the discrete second-difference operator on the positive integers. 
We standardize the Laplacian such that (A/) n := f n+ i + f n -i — 2/ n for / = (/„) € 
£ 2 (N), and consider equations of the form 

(3.1) (-A + V> = 0, 

where the potential-energy function V is a diagonal operator with real values V n . 

Eq. (|3.f [) and its solutions share many of the properties of classical Sturm- 
Liouville equations, as is laid out for example in pQ. For our purposes we recall 
that: The solution space is two-dimensional, and the Wronskian of any two solutions 

(3.2) W[u^\u^} := t#>u<2i - u£l lU W 
is constant. A Green matrix as a solution of 

(3.3) {-A + V)G = I, 

where / is the identity matrix, and every Green matrix can be written as the sum 
of a vector in the null space of (—A + V) and the particular Green matrix 

(1) (2) 

/q a\ /^(v) max(m,n) min(m,n) 

( ' m > n : " W[uW,uW] ' 

provided that {u^, u^} are linearly independent. 
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A feature of the discrete Schrodinger equation (|3.1[) that is not shared by ()2.ip 
is an invariance under the transformation 

u n -> {-l) n u n 

(3.5) V n -> -4 - V n , 

as can be easily checked. Among other things, this implies that any fact proved 
under the assumption, for example, that V n > has a counterpart for V n < —4, 
with systematic sign changes. 

Our goal in this section is to present an analogue of the Bohl transformation 
for the discrete Schrodinger equation (|3.ip . In particular, we offer a discrete version 
of some of the results of [10] . §4, and show in particular that the diagonal elements 
G nn of the Green matrix allow the full solution space to be recovered formulaically. 
We build on some earlier steps in this direction by Chernyavskaya and Sinister 
[£H \f\. As in |10j we furthermore point out connections between the diagonal of 
the Green matrix and an Agmon distance for (|3.1[) . 

In the discrete situation the use of exponentials of integrals is not the most 
natural, so we instead seek to represent a pair of solutions in the forms 

n / n 

(3-6) (fit = Zn JJ Si, ip~ = z n f JJ St 

1=1 \e=i 

Since the product of these two solutions is the diagonal of a Green matrix, up to a 
constant multiple, this suggests that if we begin by selecting a Green matrix such 
that G nn is nonvanishing, then we can directly define 




z — (G ) 1/2 

It remains to work out the most convenient form of S^ 1 . If the Wronskian is 
scaled so that W[ip~ , ip + ] = 1, then substitution of the ansatz (|3.6j) leads after a 
calculation to 



(3-7) S n - — = . 

On ZnZn—1 

Here we pause to observe two ambiguities in relating ip^ to the potential V. 
The first is that, due to the invariance fl3.5p , if 

(3'^) G mn ^niin(m,ri) ^max(m,n) 

is the Green matrix for some potential function V n , then the same diagonal elements 
G nn also belong to the Green matrix for an equation of type (|3 . 1 [) but with potential 
function V n — —4 — V n - Secondly, (|3.7[) is equivalent to a quadratic expression for 
S n , and therefore the solution is generally nonunique. These ambiguities are avoided 
when the Schrodinger operator H = — A + V in (|3.1[) is positive, so G nn > and 
by convention z n > 0. We can then fix S n as the larger root of (|3.7|) . 
Accordingly, in this situation we simply define 

1 + Jl + Azjz*^ 
2z n z n —i 



(3.9) 



0„ 
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A pair of functions (p^ can now be denned by the ansatz (|3.6|) . i.e., when expressed 
in terms of z n , 



(3.10) <?x~zn n 

k=m-\-l 



l + y / l + 4z2 z 2_ 1 
1z n z n —i 




<rr< 



Remarkably, with this definition, both (p + and tp solve a single equation of the 
form (|3. 1[> . where the potential function V n is determined from z n via 



Ap+ _ l + V 1 + 4z » z »+i 2z 2 n 



T/H — ^» _ V j ^"-1 9 

K n - — + — r, 9 ' / 

^ 2 ^ l + ,/l+4^4_ 1 



2zi 

/0 11\ _ Z n+1 „[z] | Z n-1 r, 

(^•11) A„ +1 + — ^-2, 

provided that V n > —2. (Else a different root must be chosen in (J379J) . ) To see that 
ip~ solve the same discrete Schrodinger equation, let us separately calculate 

(3.12) ^ = -^±±- + ^l5W - 2, 

and note that since Sffl has been chosen to satisfy (|3.9Jt . the difference between 
these last two expressions is 

1 1 _ 

z n z n 

This leads to a theorem in the spirit of |10j . 

Theorem 3.1. Suppose that Q3.1[) /ias two independent positive solutions for 
m < n < N, with N > M + 2, and denote the associated Green matrix G mn . 
Since G nn > for m < n < N , we may define z n := \J G nn > 0. In terms of z n , 
determine Sffl and according to (|3. 91) and (|3.10[) . TTien 

(1) (p^ is an independent pair of solutions of (|3.ip for m < n < N . 

(2) G nm = z n z m nLm+i TF\> M < m < n < N. 

(3) The potential function is determined from G nn by a nonlinear difference 
equation, 

(3.13) - (\A + 4G Tl „G„ + i„ + i 

n — 1 n — 1 

= (K + 2)G 

Remark 2. In what follows we are mainly concerned with what happens when 
N — > oo. In that case the assumption that there are two positive solutions is a 
question of disconjugacy in the theory of ordinary differential equations, cf. |20L [Tj . 
If, for example, V n > for n > No, then it is not difficult to show that no solution 
can change sign more than once, and that therefore the positivity assumption is 
satisfied for n sufficiently large. As will be seen in the proof, a necessary condition 
for the assumption is that V n > — 2. 
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Per the symmetry remarked upon in p.5[) . an alternative to positivity is the 
assumption that there are two solutions ip^ such that (—l) n ip^ > 0. A sufficient 
condition for this is that V n < —4 and a necessary condition is that V n < —2. 

Proof. The essential calculation was provided in the discussion before the 
statement of the theorem. Given that the Wronskian of ip~ and cp + is 1, these two 
functions are linearly independent and therefore a basis for the solution space of 

(-A + yWv = o, 

Vn^ being defined by (I3.11[) . Moreover, 

Gmn ^min(m,n) ^max(m,n) 
\z\ 

is a Green function for — A + Vn . 

Hence the crux is to show that Vn^ is the same as the original V n of (|3.ip . 
Because S n was defined such that 

SM _ _L - _J_ 

" eM Z Z 1 ' 
O^i ' ^n^n—1 

we may rewrite (13.111) as 

(3-14) + 2 = -L ^l + 4z ? ^ +1 + v /i + 4z^2_ i 

From the definition of z n and the assumptions of the theorem, we know that for some 
independent set of positive solutions ip^ of (|3.1[) . with Wronskian 1, = ^ip~. 
Therefore 

4^4±1 = 4(^~ ± i)(Wn±l) 

= (^+C±l+fc^±l) 2 -l- 

Hence (IXT4]l yields 



1 



V n + 2, 



as claimed, and establishes (13.13[) . □ 

It may well be asked at this stage why we have restricted ourselves to the 
situation where G nn > 0, for at the formal level the calculations given above remain 
valid without assuming positivity. In the discrete setting, continuity is not available 
to connect the values of a solution <p n as n varies, and hence without an assumption 
such as positivity, there is a degree of indeterminateness in defining solutions by a 
prescription such as (|3.9[) . For some choice of phases in (|3.14l) . it will still be true 
that Vn as defined there coincides with V n , but the implicit nature of these choices 
of square root is problematic. Possibly a suitable canonical choice of phase or ideas 
from Teschl's oscillation theory for Jacobi operators |28j could help avoid implicit 
definitions, and we hope to elaborate this point in future work. 
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Returning to the case where G„„ > 0, Formula (|3.10|) suggests that S n can be 
related to an Agmon distance O 121] . that is, a metric dA{m,n) on the positive 
integer lattice such that every £ 2 solution <fr of (|3.ip satisfies a bound of the form 

and that as a consequence <f)~ decays rapidly as n — >■ oo. Thus if z n is bounded we 
expect an Agmon distance to be something like X)™=m+i l n $l j assuming n > m. 
(We write the Agmon distance in this way because the triangle inequality is an 
equality on the integer lattice, which implies that any metric takes the form of a sum 
of quantities defined at values of I from m+ 1 to n.) In Agmon's theory, however, it 
is desirable that the distance function be a quantity that can be calculated directly 
from the potential alone (or at least dominated by some such expression). As we 
shall now see, understanding the diagonal of the Green matrix allows the derivation 
of Agmonish bounds. We begin by showing that G nn is comparable to (V n + 2) _1 
in a precise sense. 

Lemma 3.2. Suppose that liminfn^oo V n > C > 0, and let G mn be any positive 
Green matrix for (|3.1j) on the positive integers. Define 



C(C + 2)J C(C + 2)' 



f 2 \ 2 2 

Then for n sufficiently large. 
Consequently, 



y/(V n + 2)(K-1 + 2) + y/4 + (V n + 2)(K-1 + 2) 

2K A - bn 



(3 16) < V(Vn + 2)(K-1 + 2) + y^+jVn + 2)(K-1 + 2) 

Remark 3. The upper bound is of the same form as a semiclassical upper 
bound proved in |19j . To simplify it, Ka could be replaced in these inequalities by 



1 + ^>K A 



(see proof). 



Proof. The lower bound on G nn is immediate from Statement (3) of Theorem 
13.11 the left member of which is larger than 1. 

The upper bound in (|3.15[) requires a spectral estimate. The Green matrix G mn 
is the kernel of the resolvent operator of a self-adjoint realization of —A + V on 
l 2 ([N, oo)) for some N, where the boundary condition at n — N, N + 1 is that 
satisfied by ip+. Since — A > on this space (as an operator), inf sp(— A + V) > C, 
and hence, by the spectral mapping theorem, ||(— A + ^) _1 || p < C^ 1 . Since 
G nn = (e„, (—A + V)~ 1 e n ), where {e„} designate the standard unit vectors in £ 2 , 
it follows that G nn < C _1 . Inserting this into (|3.11j) would already imply (|3.15j) 
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with Ka replaced by y/l + 4/C 2 . To improve the constant, replace only the terms 
G„±i„±i in (ESP b y getting 



(3-17) (V n + 2) < v C . 

Since 

y"! + £2/ 

a; 

is a decreasing function of x when x, y > 0, an upper bound on G nn is the larger 
root of the case of equality in (|3.17l) (which is effectively a quadratic). The claimed 
upper bound with the constant Ka results by keeping one factor V n + 2 in the 
solution of the quadratic, replacing the others by C + 2. 

The bounds on sffl result from inserting the bounds on G nn into (|3.9p and collecting 
terms. □ 

We can now state some Agmonish bounds. 

COROLLARY 3.3. Suppose that liminfn^oo V n > C > and fix a positive integer 
m. Then the subdominant (i.e., eventually decreasing) solution <fi~ of (|3.1|) satisfies 
(a) 

vi + 2^ 



K A J 

(b) If, in addition, n(V n+ i — V n ) £ l l , then 



n V, + 2 + 7W4) 
\e=m ' ) 

Proof. The ansatz (|3.9p allows an identification of <\>~ with a constant multiple 
of tp~ , in the representation p.lOj) . Because z n is bounded, so is 



IK 



We then use the lower bound on S e from the lemma, but simplify by dropping the 
4, which allows the product to telescope in a pleasing way, producing (a). 

For (b) we note that the additional assumption on V n allows us to conclude that 
ip is well-approximated by a Liouville- Green expression in [19] . Theorem 4.1, which 
is a bounded quantity times the reciprocal of the expression in parentheses. □ 

Thus when liminf, woo V^ > 0, a suitable Agmon distance dA{m,n) for (|3.1[) is 
given by 

n 

(ln(VJ + 2)-ln*a), 

e= m +i 

or by 

V e + 2+ VVH^ + 4) 

provided that n(V n+ i — V n ) G I 1 . 

We close with a Darboux-type factorization for a generic discrete Schrodinger 
equation (|3.1|) . A Darboux-type factorization for general Jacobi operators was 
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previously considered by Gesztesy and Teschl in |15] . As in §2] the novel feature of 
this factorization is that it is constructed using the diagonal of the Green matrix. 

To this end, choose a Green matrix such that G nn is nonvanishing for a range 
of values of n, and define a solution <p+ according to (|3.9I) . The phase of the square 
roots is chosen (if necessary) to ensure that v}^ from (|3 . 14[) equals V n . 

Theorem 3.4. Given a Green matrix such that the diagonal Gkk is nonvanish- 
ing for n — 1 < k < n + 2, and choosing the phase of the square roots as described 
above, 



-A+V n = K 



_ l + (l + 4G ra „G»+i» +1 ) 1 / 2 ' 
2G nn 



l + (l + 4G im G n+ln+1 ) 1 /2 

where 1Z is the shift operator such that [7Zf] n = f n -i an d the right- difference oper- 
ator is defined by [V+/] n := f n+1 - f n . 



Remark 4. As with (|2.9I) . there is a second factorization, with shifts and 
differences reversed, and n + 1 replaced by n — 1. 

PROOF. Writing 

n _ 1 1 + (1 + 4G Im G„+i ri+l) 1//2 _ 1 ^n+l^n+l 

with defined in (|3.9|l . we first note that, by a simple calculation, 
(3.18) [V+ + Q n ] ip+ = 0. 

This motivates calculating 



H 



Qr, 



[V+ + Q n ] 



1-Qn 

which is well-defined because Q n ^ 1, owing to (|3.9[) with z% nonvanishing. The 
left factor was chosen to produce a convenient cancellation, ensuring that H has 
the form of a discrete Schrodinger equation, with a shifted index: 

(Hf) n = (-A/) n+1 + (y^F^ - Qn+lJ fn+1 + • /„. 

We now verify that when the index is shifted back, the potential term is indeed V n : 

Qn-l _ n \ _ ( Zn-1 _-t \ _ f-i _ Z n+l^n+l \ 
z n— 1 ^n+l^n+l 



= -2 + 

which reduces to V n according to (|3.11j) . □ 
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